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A Closed  Form  Expression  for  Line-of-Sight  Probability 


PURPOSE 


This  report  describes  a closed  form  expression  for  the  probabilitv  o^ 
line-of-sight  (Pine)  i'o»'  a given  terrain  type  as  a function  of  observer  altitude  ' 
and  observer-to-taPget  range  R. 

Figures  illustrate  the  goodness-of-fit  of  "Pine  vs  R"  curves  to  actua'' 
terrain  and  terrain  classification  type  data. 


Next  page  is  blank. 


®5^!^SI?;5'‘®t®^’^^*#?S?f*'<S''SK>Sgrr^  f**ww-<3f'>«(!‘>  x^wryyc»i>y»»‘« ^""Ky* *<<a»iMa»’^<»«-»«f nH»w«^  * 


A Closed  Form  Expression  for  Line-of-Sight  Probability 


1.  BACKGROUND 

For  several  years,  those  involved  with  modeling  of  line-of-sight  orobability 
(P|^0s^  have  been  investigating  analytic  approximations  for  specific  ard  classi- 
fication type  terrains.  The  authors  propose  the  expressions  in  this  text  as  a 
meaningful  first  step  toward  the  accomplishment  of  establishing  such  analytic 
expressions. 

2.  ASSUMPTIONS  AND  HYPOTHESES  I 

It  is  assumed  that  P,Qg  is  1.00  when  an  observer  looks  down  where  lookina  | 

down  defines  a zero  elevation  angle,  6 = 0 degrees.  The  expression  selected  | 

^ should  yield  "Pj^q^  vs.  R (observer/target  range)"  curves  that  exhibit  a slight  I 

I rounding  at  the  near  range.  That  is,  P.Qg  is  approximately  1.00  for  elevation  | 

; angles  less  than  60  degrees  (Ref.  1).  An  initial  hypothesis  is  formed  bv  taking  I 

the  exponential  of  a cotangent  (ctn)  function  on  the  elevation  angle  1 

‘ i 

5 i 

f 

Plos  = 1 - exp  (-a  ctn  0}  (1)  ’ 

t 

’ where  elevation  angle  0 is  interpreted  as  the  arctangent  of  the  ratio  of  range 

! and  observer  attitude. 

^ 6 = Arctan  (R/H)  (2) 

? The  condition  that  0 _<  P^^q^  £ 1 requires  tliat  "a"  be  positive. 

3.  DEVELOPMENT  OF  Pj^q^  EXPRESSION 

An  "a"  family  of  "Plqs  versus  6"  curves  is  shown  in  Figure  1.  The  'V 

parameter  characterizes  a rate  at  which  line-of-sight  is  lost.  This  rate 
j roughness  factor  reflects  the  combined  degradation  effects  of  terrain  curvature, 

target  height/profile,  and  vegetation. 

Comhinina  eouations  1 and  2 yield. 

P^^^  = 1 - exp  (-a  H/R)  (3) 

Figures  2 and  3 illustrate  equation  3 as  an  "H"  family  of  "Plq^  vs.  R"  curves. 

Figure  2 applies  for  a rate  roughness  value,  a = 10  while  Figure  3 applies  for 

a rate  roughness  value,  a = 20.  Maintaining  a constant  product  of  a and  fl  will 

maintain  the  same  P^^q^  value.  For  example,  if_a_is  doubled  and  II  iThali^, 

p, remains  unchanged  for  a given  range. 

Lub 
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b which  when  multiplied 
as  R approaches  infinity. 


(4) 


vs.  R"  curve,  it  is  referred 
to  as  the  asymptote  roughness  factor.  The  need'for  considering  non-unitv 
values  {0<b<l)  for  this  asymptote  factor  was  hypothesized  as  necessary  since 
real  data  appear  to  show  curves  that  approach  values  greater  than  zero 
at  far  ranges  (Ref.  3). 

Equation  4 was  fitted  to  three  specific  maps  of  the  5Dd  classification  for 
various  observer  altitudes  (Ref.  2).*  Figures  5 through  10  illustrate  that 
equation  4 fits  P.qj  data  well  for  each  specific  map.  However,  these  same 
figures  demonstrate  the  variation  in  rate  roughness  factor  g with  observer 
altitude  H.  This  variation  of  a versus  H challenges  the  hypothesis  that 
Plos  fundamentally  a function  of  elevation  angle  0.  That  is,  P^^j.  remains 
constant  for  a given  ratio  of  R and  H.  Therefore,  with  this  allowance  of  a 
false  hypothesis  on  0,  a is  fitted  wTth  H.  In  addition.  Figures  5 through  10 
illustrate  that  parameter  b is  constant  at  1.00.** 

To  characterize  the  terrain  roughness  for  the  5Dd  terrain  classification, 
three  terrain  samples  in  the  5Dd  classification  were  analyzed.*** 


Figure  4 illustrates  the  behavior  of  a factor 
by  elevation  angle  0 increases  the  Pj^q^  asymptote 
The  expression  plotted  in  Figure  4 is 

Plos  1 - exp  (-a  ctn  (b  arctan  (R/H))). 
Since  b affects  the  asymptotic  value  of  the  "P,  nr 


*This  classification  is  detailed  in  Table  1 as  extracted  from  Ref.  2.  In  a 
descriptive  sense,  5Dd  is  rough  terrain. 

**For  19  of  21  terrain  fits  (one  fit  for  each  of  seven  observer  altitudes  H for 
each  of  three  terrains),  the  ^ value  was  l.CO.  On  the  other  two  cases,  the  b 
value  was  0.99. 

***Presalection  process  of  what  terrains  were  to  be  analyzed  was  based  on  t'le 
availability  of  tapes  with  digitized  terrain  characteristics  and  the  knowledge 
that  the  terrains  selected  were  different  from  each  other  by  geographic 
location  and  general  appearance.  The  Appendix  gives  a listing  of  terrain  maos 
analyzed. 
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The  a values  were  averaged  over  these  three  terrains  at  each  observer  altitude 
H.  Then  these  a values  were  least  squares  fitted  as  a function  of  observer 
altitude  H.  Equation  5 is  the  result  of  this  fitting  process. 

a = (97.59267) 

The  resulting  Pj^q^  expression  for  the  5Dd  terrain  classification  is* 

\0S  " ®^P  (-(97.59267)  . (6) 

Figures  11  through  17  illustrate  the  fit  of  equation  6 to  the  "P|Q5  vs.  R" 
curves  of  the  three  5Dd  terrains  from  which  parameters  for  eouation  t were 
derived. 

4.  RESULTS 

The  results  from  Figure  1 are: 

a.  Larger  values  of  a imply  that  the  probability  of  line-of-sight  is 
greater. 

b.  The  rate  roughness  factor  a causes  marked  changes  in  curves  of 
^LOS  elevation  angle  when  given  values  between  0.01  and  10. 

Results  from  Figure  5 through  10 

a.  For  specific  terrain  and  observer  height,  equation  3 does  fit  data  of 
"Plos  vs.  R"  well. 

b.  Equation  3 is  not  sufficient  to  depict  the  Pj^q^.  characteristics  as  both 
a function  of  range  and  observer  height.  To  obtain  such  an  expression  involves 
fitting  the  dependence  of  the  rate  roughness  factor  a on  observer  altitude  H. 
(See  Appendix) 

a = AH®  ' 


*The  Appendix  provides  a description  of  the  programs  involved  to  conduct  a fit 
of  equation  4 to  digitized  terrain  data.  Also  provided  is  the  subseouent  analysi 
to  arrive  at  this  expression. 
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Results  from  Figures  11  through  17  are: 


a.  Variations  in  "Plqj  vs.  R"  curves  among  terrains  in  the  same 
classification  appear  large. 

b.  Equation  6 does  yield  a "P^qj  vs.  R"  curve  within  the  variation  of 
the  three  terrains  of  classification  5Dd  from  which  the  two  parameters  in 
equation  6 were  obtained. 

5.  RECOMMENDATIONS 

a.  A similar  analysis/fit  should  be  conducted  on  all  of  the  Natick 
classifications  to  form  a data  base  for  immediate  use  by  expected  value  and/or 
Monte  Carlo  simulation/wargame  models. 

b.  The  Natick  classification  scheme  could  be  improved  by  grouping  terrains 
into  categories  determined  by  a specified  variance  about  given' vs.  R" 
curves. 

6.  COMMENT 

A similar  analysis  was  conducted  on  three  terrains  of  the  Natick  Classification 
lAa  - a smooth^ table -top  terrain.  Results/fits  are  as  good  as,  if  not  better 
than,  those  shown  in  this  report  on  the  5Dd  classification.  The  corresponding 
equation  6 for  this  lAa  classification  is 


P^os  = 1 - exp  (-770.748 


(6a) 
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TABLE  I.  CLASS  INTERVALS  OF  CLASSIFICATION  DESCRIPTORS 

Maximum  Hill 

Height  (Local  Relief) 

Descriptor 

Class  Interval 

1 

0 -10  Meters 

0 

33  Feet 

2 

10  -30 

33 

- 

90 

3 

30  -50 

99 

• 

165 

4 

50  -100 

165 

- 

330 

5 

100  -300 

330 

• 

990 

6 

OVER  300 

OVER 

990 

Modal  Hill 

Height  (Local  Relief) 

Descriptor 

Class  Interval 

A 

0 -10  Meters 

0 

33  Feet 

B 

10  -20 

33 

- 

66 

C 

20  -35 

66 

• 

115 

D 

35  -50 

115 

• 

165 

E 

50  -75 

165 

- 

248 

F 

75  -100 

248 

• 

330 

G 

100  -125 

330 

- 

413 

H 

125  -150 

413 

- 

495 

I 

150  -175 

495 

• 

578 

J 

175  -200 

578 

- 

660 

K 

OVER  200 

OVER 

660 

Number  of  Pcsitive  Features  Per  Mile 

Descriptor 

Class  Interval 

a 

0 -0.8/Kilometer 

0 

0.5/Mil( 

b 

0.8  -1.6 

0.5 

- 

l.C, 

c 

1.6  -2.4 

1.0 

- 

1.5 

d 

2.4  -3.2 

1.5 

- 

2 0 

e 

3,2  -4.0 

2.0 

2.5 

f 

OVER  4.0 

OVER 

? . 5 

Next  page  is  blank. 


13 


APPENDIX 


A.l.  PROCESS  TO  OBTAIN  EXPRESSION 


Obtaining  a closed  form  expression  for  Pj^q^  involves  three  stens: 

Step  1:  Generating  P.q^  curves  from  digitized  terrain  da^a  and  fitting 
equation  4 to  these  curves. 

Step  2;  Obtaining  a correlation  function  between  rate  roughness  factor 
a and  observer  altitude  H,  and 

Step  3:  Calculating/plotting  of  the  theoretical/fitted  curve. 

Step  1 consists  of  the  program  "real  curve  fit."  This  program  is  a 
combination  of  two  programs.  The  primary  program  "LOS  Rings,"  written  bv 
Warren  Olsen,  generates  the  "^LOS  vs.  R"  curves  from  digitized  terrain  data 
(Ref.  3).*  The  secondary  program  "Minmax"  modeled  by  J.  Sheldon  and  written/ 
executed  by  Sean  Smith,  fits  equation  4 to  the  digitally  processed  "P,^-  vs.  P" 
curve  by  using  the  Kolmogorov-Smirnov,  maximum  difference  test  for  Goodness-of- 
Fit. 

The  program  "Real  Curve  Fit"  is  written  in  FORTRAN  IV  compatible  only 
with  the  BRLESC  II  computer  at  Aberdeen  Proving  Ground,  Maryland.  Table  II 
illustrates  the  information  necessary  to  run  a single  case  - a case  consists  of 
a "Plq5  vs.  R"  fit  for  each  of  seven  observer  altitudes. 

Table  II 

Computational  Information  on  BRLESC  II  Program  - Real  Curve  '^it 


Run  Time 
(Min) 

Range 

4 G 8 

45 

71 

105 

Output 

(Lines) 

600 

735 

870 

*The  P|Q^  result  from  digitized  terrain  is  an  average  result.  The  averagino 
occurs^over  observer  location  on  the  terrain  over  a 360  degree  scan.  This 
process  accounts  for  the  somewhat  smooth  curve  representing  th<»  actual 
"Plos  vs.  R"  result. 
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step  2,  finding  an  equation  which  relates  rate  roughness  factor 
a to  observer  altitude  H,  uses  a standard/library  program  "Regression 
Analysis  with  Plotting.^  This  process  begins  by  averaging  all  the  rate 
roughness  factors  a over  the  three  terrains  for  each  of  the  observer 
heights  H.  These  average  values  of  a are  least  squares  fitted  to  H 
in  the  form* 

a = AH®.  (7) 

Step  3 consists  of  the  plotting  of  equation  6 which  is  the  result 
of  combining  equations  3 and  5.  The  general  form  for  equation  6 is 

P^os  = 1-  exp  (-AH^®^^Vr)  (8) 

A. 2.  TABLE  OF  ACTUAL  TERRAINS  STUDIED 

Table  III  lists  the  terrains  studied  in  this  report.  These  terrains 
are  shown  in  Figure  A.l. 


TABLE  III 

Actual  Terrains  Studied 


Natick  Terrain 
Classification 

DMA  m 
Numbers 

Actual  Terrain  Code**/City 

5Dd 

L6134  - L6132 

G-93/Bayreuth 

L5120  - L5118 

G-4C/Stadt  Allendorf 

L4920  - L4918 

R-44/Bad  Wildungen 

lAa 

L7130  - L7128 

M 

G-70/Nordl ingen 

L7332  - L7330 

G-76/Neuberg 

L7334  - L7332 

G-77/Ingolstadt 

\ 


*The  correlation  function  between  4.  and  J1  was  tested  in  two  different  forms 
(a=AH“  and  a=AB").  The  best  results  were  obtained  by  using  the  former 
expression,  Equation  7. 

♦♦Terrain  codes  are  standard  DMA  (Defense  Mapping  Agency)  codes  (Ref.  4). 
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Figure  4.  A Family  of  P^Qg  v; 

1-exp  (-0  ctn  (b  artan  (R/H)  ) ) 


OBSERVER  HEIGHT  »3.00M 

TARGET  HEIGHT  = .00 

AREA  G93 

CLASSIFICATION  5Dd 

RATE  SHAPE  FAQOR  a * 43.11 

ASYAAPTOTIC  SHAPE  FAQOR  B * 1.00 

MINIMUM  AAAXIMUM  DIFFERENCE  « 0.03 


Figure  5 Une-of- Sight  Probability  vs.  Range. 

(Artificial  Girve  Fit) 


i OBSERVER  HEIGHT  » 1.67  M 

TARGET  HEIGHT*  ,00 
j AREA  G93 

1 CLASSIFICATION  5Dd 

1 RATE  SHAPE  FACOR  a » 61.96 

ASYMPTOTIC  SHAPE  FACTOR  B « 1.00 
! MINIMUM  AAAXIMUM  DIFFERENa  > 0.02 


Figure  e Une-of- Sight  Probobiiity  vs.  Range. 

(Artificial  Grve  Rt ) 


RANGE  (KM) 

OBSERVER  HEIGHT  » lO.OOM 

TARGET  HEIGHT*. 00 

AREA  G44 

CLASSIFICATION  50  d 

RATE  SHAPE  FAaOR  a , 40.76 

ASYMPTOTIC  SHAPE  FAaOR  B * 1.00 

MINIMUM  AAAXiMUM  DIFFERENCE  * 0.05 


Figure  7 Une-o^- Sight  Probability  vs.  Range. 

(Artificial  Gjrve  Fit ) 


OBSERVER  HEIGHT  » lOO.OOM 

TARGET  HEIGHT*. 00 

AREA  G48 

CLASSIFICATION  SOd 

RATE  SHAPE  FACTOR  a * 2767 

ASYMPTOTIC  SHAPE  FACTOR  B » 1.00 

MINIMUM  AAAXIMUM  DIFFERENCE  * 0.03 


Figure  8 Line -of- Sight  Probability  vs.  Range. 

(Artificial  Grve  Fit ) 
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RANGE  (KM) 

OBSERVER  HEIGHT  » 500.00M 

TARGET  HEIGHTS. 00 

AREA  G93 

CLASSIFICATION  5Dd 

RATE  SHAPE  FAaOR  a*  15.36 

ASVmOTIC  SHAPE  FAaOR  B * 0.99 

MINIMUM  AAAXIMUM  DIFFERENCE  * 0.03 


Figure  9 Line -of- Sight  Probability  vs.  Range. 

(Artificial  G^rve  Fit ) 
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RANGE  (KM) 


OBSERVER  HEIGHT  « 20.00  M 

TARGET  HEIGHT*. 00 

AREA  G44 

CLASSIFICATION  50  d 

RATE  SHAPE  FACTOR  a » 29.70 

ASYAAPTOTIC  SHAPE  FACTOR  B * 1.00 

MINIMUM  AAAXIMUM  DIFFERENCE  « 0.06 


Figure  io  Line -of- Sight  Probability  vs.  Range. 

(Artificial  Grve  Fit ) 


<SSteuW!h>9lf<SfW  '•VSH\l!*iiV9^<'>ifvyMuw*¥*t}rtyeM^  >^W?^^TS^!??B8Sg^8S8SSBySgg^?W^1^y^^^  iff®l^^^^?13®a?3 


RANGE  (KM) 


OBSERVER  HEIGHT  » 1.67M 
TARGET  HEIGHT  = .00 
ROUGHNESS  FACTOR  a « 83.24 
ROUGHNESS  FACTOR  B = ].00 
OASSIFICATION  5Dd 


Figure  ]]  Line -of -Sight  Probability  vs.  Range. 

(Artificial  Curve  Fit) 
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RANGE  (KM) 


OBSERVER  HEIGHT  * 3.00  M 
TARGET  HEIGHT  « .00 
ROUGHNESS  FACTOR  a « 69.39 
ROUGHNESS  FACTOR  B»1.00 
CLASSIFICATION  5D<1 
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OBSERVER  HEIGHT  * 10.00  M 
TARGET  HEIGHT  * .00 
ROUGHNESS  FACTOR  a » 47.75 
ROUGHNESS  FACTOR  BM.OO 
OASSIFICATION  5Dd 


Figure  13  Line -of -Sigh I*  Probability  vs.  Range. 

(Artificial  Curve  Fit) 
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i 0 1 2 3 4 5 

I 

RANGE  (KM) 

1 

j OBSERVER  HEIGHT  * 20.00  M 
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Figure  14  Line-of- Sight  Probability  vs.  Range. 

(Artificial  Curve  Fit) 
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OBSERVER  HEIGHT  = 50.00  M 
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CLASSIFICATION  50  d 


Figure  Line-of- Sight  Probability  vs.  Range. 

(Artificial  Curve  Fit) 
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OBSERVER  HEIGHT  » 100.00  M 
TARGET  HEIGHT  * .00 
ROUGHNESS  FACTOR  a * 23.37 
ROUGHNESS  FACTOR  B»1.00 
CLASSIFICATION  50  d 


Figure  Line-of- Sight  Probability  vs.  Range. 

(Artificial  Curve  Fit) 
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OBSERVER  HEIGHT  * 500.00  M 
TARGET  HEIGHT  * .00 
ROUGHNESS  FACTOR  a * 14.18 
ROUGHNESS  FAaOR  B = 1.00 
OASSIFICATION  50  d 


Figure  Line-of- Sight  Probability  vs.  Range. 

(Artificial  Curve  Fit) 
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Figure  A-1.  LocoHon  of  Tenroir)  Areas. 


